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Using non-equilibrium Green functions we calculate the spin-polarized current and shot noise in 
a ferromagnet-quantum-dot-ferromagnet (FM-QD-FM) system. Both parallel (P) and antiparallcl 
(AP) magnetic configurations are considered. Coulomb interaction and coherent spin-flip (similar 
to a transverse magnetic field) are taken into account within the dot. We find that the interplay 
between Coulomb interaction and spin accumulation in the dot can result in a bias-dependent current 
polarization p. In particular, p can be suppressed in the P alignment and enhanced in the AP case 
depending on the bias voltage. The coherent spin-flip can also result in a switch of the current 
polarization from the emitter to the collector lead. Interestingly, for a particular set of parameters 
it is possible to have a polarized current in the collector and an unpolarized current in the emitter 
lead. We also found a suppression of the Fano factor to values well below 0.5. 

PACS numbers: PACS number 



I. INTRODUCTION 

Spin-dcpcndcnt transport in quantum dots is a subject 
of intense study nowdays due to its relevance to the new 
generation of proposed spintronic devices that encom- 
passes, for instance, the Datta-Das transistor^ memory 
devices^ and as an ultimate goal quantum computers.— 
In particular, the recent progress in the coherent con- 
trol of electron spins in quantum dots^i 7 - has stimulated 
even further the research in this field, for possible appli- 
cations in quantum computation and quantum informa- 
tion processing^ In addition to these fascinating tech- 
nological applications, quantum dots constitute a unique 
well-controllable system to study fundamental physical 
aspects of transport in the strong Coulomb-correlated 
regime, and its interplay with spin-dependent effects. 

A common geometry used for transport studies in 
quantum dots consists of two leads weakly coupled to 
a quantum dot (QD) via tunneling barriers. Spin- 
dependent effects such as spin accumulation and spin- 
polarized transport can occur in these systems when 
both leads are (or at least one of them is) ferromagnetic 
(FM). The junction FM-QD-FM resembles the standard 
TMR&i£ and GMRii geometries composed of an in- 
sulator layer sandwiched by two ferromagnetic metal- 
lic leads, except for the quantum dot replacing the in- 
sulator layer. This system (dot coupled to FM leads) 
was recently experimentally realized in the context of 
semiconductor quantum dot a 12 i 13 and molecules i 14 ' 15 i 16 
A wealth of novel spin-dependent effects has been ob- 
served in this system due to the interplay of quan- 
tum confinement, Coulomb correlations, Pauli principle 
and lead-polarization alignments. For instance, novel 



effects such as spin- accumulation , 7 i 8 spin-diode , ° 
spin-blockade ) 21 ' 22 i 23 ' 24 spin-current ringing ) 25 i 26 nega- 
tive differential conductance and negative TMR "' 21 and 
so on, arise in this context. In order to obtain additional 
information, not contained in the average current, shot 
noise has also been analyzed in several spintronic sys- 
tems. A few exemples include shot noise in spin- valve 
junction o 27 ' 28 i 29 ' 30 and quantum dots attached to ferro- 
magnetic leads i 31 i 32 i 33 

Here we apply the Keldysh non-equilibrium technique 
to study spin-polarized transport (current and shot noise) 
in a FM-QD-FM system (Fig. 1). Both parallel (P) and 
antiparallcl (AP) lead magnetization alignments are con- 
sidered. The left and the right lead materials are taken to 
be different, thus resulting in additional effects, not seen 
for leads with the same material. We analyze both the 
current and the shot noise in the presence of Coulomb 
interaction and spin- flip in the dot. We find an inter- 
play between spin accumulation and Coulomb interaction 
that gives rise to a bias dependent current polarization 
p. More specifically, p can be suppressed or enhanced 
and have its sign changed, depending on the magnetic 
alignment and the bias voltage. We also note that the 
spin-flip can switch the current polarization as it flows 
from the emitter to the collector lead. In particular, it 
is possible to have an unpolarized emitter-current and a 
polarized collector-current. For the shot noise, we find 
that spin-flip can suppress it (in the AP case) with Fano 
factors reaching well below 1/2. 

The outline of our paper is as follows. In Sec. II we 
describe in detail our model Hamiltonian. In Sec. Ill we 
present the current and the noise calculations, respec- 
tively, including general formulas for these quantities. In 
Sec. IV we present and discuss numerical results for the 
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FIG. 1: Schematic of the system studied: a quantum dot 
coupled to two ferromagnetic leads. On the forward bias the 
electrons tunnel from the left lead to the right lead via the 
quantum dot. We consider configurations in which the FM 
leads are parallel (P) or antiparallel (AP). The dot has a sin- 
gle orbital level and can hold at most two electrons of opposite 
spins. Intradot Coulomb interaction and spin-flip are consid- 
ered. 

current and the shot noise. We summarize our conclu- 
sions in Sec. V. Technical details of our calculation are 
described in the appendixes A-D. 



II. MODEL SYSTEM AND HAMILTONIAN 

Our system consists of a quantum dot with one quan- 
tized level coupled to two ferromagnetic leads via tun- 
neling barriers. While the electrons in the leads are non- 
interacting, the electrons in the dot experience Coulomb 
repulsion and spin-flip scattering. The system Hamilto- 
nian is 

H = H L + H R + H D + H T . (1) 

The first three terms in ([I]) correspond to the three differ- 
ent regions: left lead, right lead and dot. The last term 
Ht hybridizes these three regions thus allowing electrons 
to tunnel from one region to the other. This term gives 
rise to current in the presence of a bias voltage. 

More explicitly, we have for the ferromagnetic leads 

H v = ek<TT]Ck<Tn C kcrTi ; (2) 

fc(7 

where eka n = e/,17 + (— 1) 1 A (Stoner Model) is the spin- 
dependent energy of the electron in lead rj = (L,R), 
with the band spin splitting A, a =|, j and #pmn = 1, 

^iT(Ti) = 0- ^he operator Ckarj i c \ arj ) destroys (creates) 
an electron with wave vector k and spin a in lead rj. The 
dot Hamiltonian is 

H D = Y^ e ° d l d <? + Un ] -n l + R(d\d l + djrf T ), (3) 

(7 



where e CT is the dot level and d a (d' a ) annihilates (cre- 
ates) an electron in the dot with spin a. Our model 
assumes a single spin-degenerate orbital level in the dot, 
e| = ej_ = e. More specifically, our dot can be singly 
occupied by an electron with spin up or down or dou- 
bly occupied by two electrons with opposite spins. We 
account for the Coulomb interaction in the dot via the 
Hubbard term with correlation parameter U. We assume 
a linear voltage drop across the system: e = eo — 
where e > 0, V is the applied voltage and e is the dot 
level for V = 0. The left \xl and the right chemical 
potentials are related by jiL — jJ-R = eV. Here we assume 
that hl is constant and defines the origin of the energy. 
For positive bias > Hr) the left lead is the electron 
emitter and the right lead is the collector. The last term 
in accounts for a coherent spin- flip in the doti^i This 
term can represent, e.g., a local transverse magnetic field 
that coherently rotates the electron spin, which can be 
experimentaly realized via ESR techniques^ or by the 
Hanle eSect^L 

Instead of carrying out the calculations with the 
Hamiltonian we perform the following canonical 
transformation, 

V 1 1=1,2 

With (|2|), the dot Hamiltonian becomes 

Hd= J2 [ £ ° + + Unm 2 , (5) 

i=l,2 

where rij = d\di. Note that in © the dot level is split into 
two levels; t\ = £q — R and £2 = £o + ^- We note that this 
canonical transformation rotates the spin quantization 
axis (e.g., to the direction of a local transverse magnetic 
field) , thus replacing the spin-flip term by a diagonal term 
with a split level. 

The tunneling Hamiltonian in |T]) is 

h t = J2{tk*ci an d a + t* ka dlc kav ), (6) 

where the matrix element tka- connects an electronic state 
in lead rj to one in the dot. Observe that the hopping pro- 
cess between the leads and the dot is spin conserving, i.e., 
tka- does not mix different spin components. Applying the 
transformation @ in ^ we find 

(— l)^"i 

ht = Y1 — -7^— {*L4 c k<"j + w c L, d 0- (7) 

karji 

Next we calculate current and noise for the model de- 
scribed above. 

III. CURRENT AND NOISE 

The current is calculated in the standard way from 
the definition Itif) = ( r L(t)) where I Lit) = —cNl is 
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the current operator, with Nl — J2k<? c k<rL c k<rL being 
the total number operator, and (...) is a thermodynamic 
average. From the Heisenberg equation Nl = JVl], 
we find 

I L {t) = ieJ^ttocL, (*)£*,(*) - C4(*)W*)L (8) 



where the labels r, a and < mean retarded, greater and 
lesser, respectively. The calculation of the retarded G\-, 
and lesser Gfj dot Green functions is presented in the 
Appendix A. 

Using this result in (|10[) we arrive at 



which results in the following current expressio: 



,38.39 



2eRe^ J 



I L {t) =2eRe^t fc0 

ka 



(t)d a (t)). 



(9) 



where 

„«,«) 



A similar expression holds for the right lead current 
Ir = —e(Nji). Since we are in stationary regime we have 
simply I i, = —Ir. Using the canonical transformation in 
Eq. ([9]) we obtain 



I L = 2eRc^t fe( 



(-1)' 



(10) 



kai 



Where <^Ll(M) 



J kL 



„{t)di(t)) is the lesser Green 
function, which is calculated via the Kcldysh nonequilib- 
rium technique! 40 ' 41 

As a starting point we construct the complex time 
Green function G ikaL (T,T') = -i(T c d i (r)4 crL (r')), 
where T c is the contour time-ordering operator and r 
and t' are complex times running along a complex 
contour. 40 ! 41 Then we go from the Heisenberg to the in- 
teraction picture, by introducing the 5-matrix operator 
S = e - l ^ dT " T{T \ Here the tilde means that H T is in 
the interaction picture. After expanding S we find^ 



x {G^(i,ti)S^.(t 1) t) + G<.(t,i 1 )S^(i ll i)} (13) 

^ (< ' o) (ii,t)=|E fe J^| 2 (-l) (i+J ' )a ^ 
9kaL '' > with the lesser Green function 
fffeo-i^i'*) = *(cj £[J L( i )cfe ( Ti( t i))' and tnc advanced 

0ne ffLi (*!'*) =^(*-*l)({ffc<ri(*l)j4<7L (*)})' hcre the 

curly brackets denote an anticommutator. 

In the steady state regime the Fourier transforms of 
the Green functions result in single frequency Green func- 
tions. Since this is the regime of interest here, we state 
for later use the Fourier transform of the leads Green 
functions. 



9kcrrj( e ) ~ 9kari( e ) 



1 



(14) 



GikaL (T , T ) = t%„ -j= 



e - eka7 1 - iS 
9k* v (t) = 2mnr,{e)8{e-e k(Tn ), (15) 

where is the Fermi distribution function of lead 77. 

A. Average current in the stationary regime 

In a stationary regime all of the Green functions de- 
pend on only t — t±, yielding the Fourier transform 

I L = 2eRe| g ]p<%(e)Ej;<(e) + G£(e)Ejf (e)}, 



x /^(r.nWny), (ID = ie / | r &{T L [(G--G> £ + G<]} l 

J c 



where Gy (r, n ) 



i(T c d t (T)d](Ti)) and gkaLin.r 1 ) 



with 



-i(T c CkaL(ri)cl a . L (' r ')}- Note tbat while G^(t,ti) is in 
the Heisenberg picture, gkah{Ti,T') is in the interaction 
picture (denoted by the tilded operators). This "sepa- 
rability" of the interaction and Heisenberg pictures fol- 
lows from the assumption of noninteracting electrons in 
the leads. This allows us to put the "difficult" part of 
the analysis entirely in the dot Green functions, which 
contain the Coulomb interaction, the spin-flip and the 
coupling to leads. 

The next step is to apply Langreth's analytical contin- 
uation rules 4 ^ to (fTTj) . to find the lesser Green function 
appearing in (jTTJJ) . This yields 



1 / + - r£ 



(16) 



(17) 



where = 27r^ fc |ifc CT | 2 £(e — tkah) is the lincwidth func- 
tion. In what follows we neglect the energy-dependence 
of F^ (wideband limit), which will be taken as a constant 
phcnomcnological parameter . 



B. Spin- resolved currents 



From Eqs. ([16]) and (|17|) we can also determine the 
spin-resolved components of the average current 



dti x 



2tt 1 2 



1 



(-If 



{GUt, h)g< aL (h,t') + G< (t, h^ihJMU) 



x [n L (G r -G a ) + G<]}. 



(18) 
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A similar result holds for If,. Equation ([18]) gives the 
spin-polarized current components with their polariza- 
tion axes defined along the magnetic moment of the leads. 
In the present study no spin-torque is considered, which 
makes the projected current the relevant quantity to in- 
vestigate. In the presence of spin-torque more general 
definitions for spin-rcsolvcd charge currents and spin- 
currents should be used. A general expression for the 
spin-current in the presence of spin transfer was recently 
derived in Ref. (H. 



C. Noise definition 

Fluctuations of the current are interesting because they 
can give additional information about the system beyond 
that provided by the average current alone.— Here we 
derive an expression for the current fluctuations, which 
include both thermal and shot noise. The thermal noise 
is related to fluctuations in the occupations of the leads 
due to thermal excitation, and it vanishes at zero tem- 
perature. Shot noise is an unavoidable temporal fluctu- 
ation of the current due to the granularity of the elec- 
tron charge. It is nonzero only for finite bias, i.e., it 
is a nonequilibrium property. In the linear response 
regime the fluctuation-dissipation theorem holds, yield- 
ing the relation S(ui) = 4fcsTG(w), where G(oj) is the 
conductance^ Hence, in equilibrium the noise contains 
the same information as the conductance. Away from 
equilibrium this relation is no longer valid and the noise 
spectrum can provide additional information. 

We define noise via S v ,y = ({SI n (t), <JJ,,/(t')}), where 
SI v (t) = Ir)(t) — I v is the current fluctuation at a time t 
in lead 77. Equivalently, 



tions, 



S w (t,t') = ({i n (t)J n ,(t')})-2ll 



(19) 



where we use the fact that = (I v (t)) = (I v i(t')) in the 
stationary regime. Using the current operator I v [Eq. 
©] and Eq. gj in Eq. JT^, we obtain 



S nn >(t,t') = (ie) 



f-' 2 

kk' (7(7' ij 



{^^ CT '(cL 7) (*)^(*)cL' 7) '(* , ) d i( i 0) 

W (4 (t)c kav (t)c{ l(r , n , ffldj (*')) 
+tt a t%, a , (4 {t)c kar , {t)d] (*')c fcW (f )>} 
+h.c. - 2/2. (20) 

D. Noise in terms of Green functions 

Each (...) term in (f2T)]) can be expressed in terms of 
a Green function. Defining the two particle Green func- 



g[ 2 \r,r') = l 2 (Tj kan (r)d l (r)c{ lalv ,(r')d J (r^ 
gi 2) (T,r') = l \Tj ka7l (T)d l {T)d){T')c k , a , 11 ,(T>)) 



gi 2 \r,r') = f ! (T c dJ(r) Cfe ^(r)4, ^VW)) 
gf\r,r') = * 2 (T4(r) CfcCT? ,(r)4(r') Cfc ' CT '^(r')), 

we can write ([20)) as 

S w (t,t>) = e 2 \{-l) iS ^{-l) jS ^ 

kk' aa' ij 

{t k at k/a ,g[ 2)> (t,t') 
-t krr t* k , a ,gi 2)> (t,t') 

-*LWff$ a)> (*,0 

+tt a t* k , a ,g {2)> (t,t')} 

+h.c. - 2/, 2 , (21) 

(2) > 

where g\ (t, t') is obtained from the complex time 

Green function g\~ (t, t') via analytical continuation. 
Similarly to the current calculation where we de- 
velop an 5-matrix expansions in G aka h(j-,T IS ) to obtain 

G^ k(xL (t,t') , here we expand the S'-matrix in (t, t') 

( 2) > 

and then obtain g\ ' (t,t'). This procedure follows the 
standard calculations proposed in Ref. [39[ to derive the 
current equation. The details of this S'-matrix expansion 
are presented in Appendix B; here we simply state the 
results, 

9i\?y) = \ E (-ir^-ir 5 ^ 



i-, i 2 — 1,2 



xti„t 



k<j b k' 



dTidT 2 gka V (n,T)g k 'a' V '(T 2 ,T') 



{G lil (r, ri)Gj l2 (t', t 2 ) - G ll2 (t, T 2 )G jil (r', n)}, 



(22) 



and 



ff2 2) ( T , T ') = -Sk<rri,k'cr'r)>gk<TT,(T' ,T)Gij(T, t') 

+ \ E (-i) <i ^(-i) <a4 " 1 

ii ,?2— 1,2 

x*L*fcv J J dT 1 dT 2 g kari {T 1 ,T)g k , (r' ,r 2 ) 

{G in (r, Ti)Gi 1J (r 2 , r') - Gy(r, r')Gi 2il (r 2 , n)}. 

(23) 

Equations (|22"|) and (|23"1) hold on a Hartree-Fock or other 

mean- field theory (sec details in Appendix B). The other 

(2) (2) 

two Green functions g 3 and g\ are given by 



ff3 (T.' r ) = 92 (t,t) 



(2)* 



(24) 
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and 



(2)/ /\ (2)*/ / \ 

9i (r,T) = g\ (t,t ). 



(25) 



From a diagrammatic point of view the terms in Eqs. 
([22]) and ((23]) involving 

Gin (r, n)gkan(n,T)Gj i2 (r', T 2 )gk'a' V ' fa, r'), (26) 
and 

G aa {r, n)gkar,(ri,T)g k 'cr' v '(T' , T 2 )G a > a > (r 2 , r') (27) 
arc disconnected. These disconnected terms, together 

(2) 

with similar ones in the equations for (r, r ) and 
<74 (t, t'), cancel identically the term 2/ 2 in Eq. (|2T|) 
(see Appendix C). So we can say that this corresponds 
to the linked cluster expansion to the noise. The other 
terms in Eqs. (|2"2"|) - (f23]) give the connected diagrams and 
thus can give a contribution to the noise. Substituting 
the connected terms of Eqs.([22] ) -(|25 |) in Eq.(|2T]) we find 



5 w (t,0-e 2 ^|t fc<T | 2 {^^(-i; 



x[g> aT] (t,t')G<(t',t) + G>(t,t>)g< av (t',t)] 
-e 2 M'lWl 2 / f d Tl dT 2 

1.1.' J c J c 

fifi (7(7 

X {Gu 2 (t, Ti)gu'a'r)' (t"2, t )Gji 1 (t , T\ )gka n (t"i, t) 

—Gij(t, t)gk>a<n'{t ,T 2 )Gi 2 i 1 (t 2 , n ) g fe(TI) (ri , i) 

— gka n {t,Tl)G ill2 (ti, T 2 )gk'(7' n '{T2 1 t)Gji(t , t) 
+gkar,(t,T 1 )G li:j (T U t')g k ( a , rr (t',t 2 )G l2l (T 2 ,t)} t>t ' 

+ h.c, (28) 

where the superscript t > t' means that an analytical 
continuation should be performed by applying Langreth's 
rules. 

E. Zero-frequency shot noise 

The shot noise is defined as the Fourier transform of 
Sr)ri'{t — t'), which in the stationary regime reads 

/>oo 

S w (w)= d(t-l/)e i ^ t -^S w (t-t'). (29) 



Using the analytical continuation of Eq. (|28|) into Eq. 
(f2"9"]l we find the following zero- frequency shot noise ^ 



e 2 r 

Sll{oj = 0) = — / deTr{m L r i G > 



- i(l - n L )r L G< +r L G > r i G < 

+ r L (G r - G a )T L (n L G> - (1 -n L )G<) 

- n L (l - n L )(G a T L G a T L + G r T L G r T L )}, (30) 



where G> satisfies the identity G> G< = G r - G a . 
All the Green functions in Eq. ([3D]) are in the frequency 
domain. In our analysis we take only the component r\ = 
jj = L. Since the dc noise is position independent, we 
have simply S LL (0) = S RR (0) = -S LR (0) = -S RL (0). 
Equation (|30|) can be expressed in a standard form as 
follow^^ (Appendix D), 

Sll(u = 0) = ^ I de 

x Tr{[n L (l - n L ) + n R (l - n R )]T{e) 



[riL 



■nH) 2 T(e)[l-T(e)]}, (31) 



with the transmission matrix T = r L GT fl G a . In the 
calculation leading to Eqs. ([30| or (f3Tj) we have truncated 
an S-matrix expansion by breaking two-particle Green 
function into products of one-particle Green functions. 
This procedure holds in a mean field theory. Thus for 
a consistent application of Eqs. (|30[) - (j3"T]) a similar ap- 
proximation (Hartrcc-Fock like) for the Green functions 
should be made (see Appendix A). Some limitations im- 
posed by this approximation are discussed in the end of 

sec. nm 



F. Model for the FM leads 

The fcrromagnctism of the leads is considered via the 
spin-dependent parameter T^. From the Stoner model, 
for instance, we can see that the density of states for 
spin up electrons of the lead is shifted with respect to 
that of the spin down electrons. Since TJJ contains in- 
formation about the spin-dependent density of states it 
is expected that T*\ ^ Following Ref. [H we de- 

fine = r [l + (—l) S "i-pL]- The parameter T gives 
the strength of the lead-to-dot coupling and ph is a pa- 
rameter describing the degree of spin polarization of the 
left lead^ Note that Tf > Tf for p L ^ 0. This means 
that the population for spin up around the Fermi energy 
in the left lead is greater than the population for spin 
down. Similarly, for the right ferromagnetic lead we as- 
sume = r [l + (-l) s °ip R ] for the parallel (P) and 
= r [l - (-l) s ^p R ] for the antiparallel (AP) lead 
alignments. Note that for the P case we have = 
and for the AP configuration = T§, with a being the 
opposite of <t. In the present work we mostly discuss the 
Pl 7^ Pr case, i.e., a geometry in which the left and right 
leads are composed of different materials (Ni and Co, for 
instance). 



G. Numerical procedure 

The numerical results are obtained following a self- 
consistent procedure. We calculate the average 

(hij) = (4dj) = I ^G<(e), (32) 
J Zni J 
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FIG. 2: Spin-resolved occupations n-f and nj. and the spin 
accumulation, m = n-\ — nj,, against bias voltage in both 
parallel and antiparallel FM lead alignments. In the P case 
n-f < ri| while in the AP alignment n-f > nj.. These im- 
balances give rise to spin-dependent population suppressions 
in the bias range 1-3 meV. This translates into a tunable m 
with the bias voltage. Parameters: pl = 0.23, pr = 0.35, 
k B T = 200 peV, T = 10 fieV, e = 0.5 meV and U = 1 meV. 



sclf-consistcntly with Eqs. (|A9|) - (|A15|) . When converged 
solutions for the expectation values {d\dj) and dot Green 
functions are found, we determine the current [Eq. (|18[) ] 
and the noise [Eq. (|5D|) ]. This iterative schema is per- 
formed for each bias voltage. 



IV. RESULTS 
A. Spin resolved electronic occupations 

Occupations. Figures [5]Ja)-(b) show the spin up and 
spin down occupations of the dot for both parallel (P) 
and antiparallel (AP) alignments. In the P case the dot 
has a net spin-down polarization with n-f < nj., while 
in the AP configuration nt > nj_. These spin imbal- 
ances in the dot can be easily understood in terms of the 
tunneling rates TJJ adopted. The parameters p^ = 0.23 
and pr = 0.35 used here 5 -^ yield the following tunnel- 
ing rates in the parallel case: Ti = 12.3 /xeV, rf = 7.7 
fieV, Tf = 13.5 ^eV and rf = 6.5 ^cV. In the AP case 
the values of the tunneling rates to the right lead are 
swapped (F R <=^ rf ). From these rates we conclude that 
in the P case a spin up electron leaves the dot faster 
than it comes in. The opposite happens for a spin down 
electron. The imbalance of these in/out tunneling rates 
results in a larger spin down occupation in the parallel 
case, i.e., n-f < nj., Fig. Ufa). By the same token, in the 
AP alignment we have n-f > nj, as seen in Fig. E^b). 

Spin accumulation. In Fig. (Hc)-(d) we show the spin 
accumulation (m = nj — tojJ as a function of the bias volt- 



age. In the zero bias limit m is essentially zero. When 
the bias increases the spin accumulation in the P case 
assume negative values. In contrast, in the AP align- 
ment m is enhanced. In particular, in the bias range 
corresponding to a singly-occupied dot (1-3 meV)£2 the 
additional suppression [Fig. [DJc)] or the enhancement 
[Fig. [21(d)] of m is due to the spin- dependent population 
suppression that takes place in the presence of Coulomb 
interaction and spin accumulation. More specifically, in 
the AP case due to Coulomb interaction n-f tends to sup- 
press more strongly n^ than otherwise. This translates 
into an enhancement of m. In the P alignment the spin 
up occupation n-f is more suppressed than , thus m be- 
comes more negative [Fig. [U(c)]. We emphasize that this 
effect happens for both the P and AP alignments because 
we assume pl ^ Pr- For equal leads we find n-f = in 
the P case so that m remains zero in this configuration. 



B. Current and its polarization 

Current. Figures [31(a)- (b) show the current in the P 
and AP cases for U = 1 meV and R — 0. Similarly to 
the occupations, some features of the spin- up and spin 
down currents can be understood in terms of the tun- 
neling rates. For instance, their saturation values (sec- 
ond plateau) can be easily calculated from the standard 
expresssior 



,53 



r^r 



(33) 



which gives H > 1+ and H < I* in the P and AP cases, 
respectively. For the first plateau Eq. ([55)) is not valid 
and these inequalities can changed 

In the P case [Fig. H[c)] if is more strongly suppressed 
than due to the interplay of spin accumulation (n-f < 
ri|) and Coulomb interaction. This results in a supression 
of the current polarization [p = (I^—I^)/ (If + 7j_)] in the 
range 1-3 meV. On the other hand, in the AP case [Fig. 
[3](d)] 1^ is more suppressed than ij. due to the inverted 
inequality nj > n^, thus resulting in an enhancement of 

P- 

Spin-flip effects. In figure ^ we show the current po- 
larization against bias voltage for distinct spin-flip pa- 
rameter R. The polarization is calculated for the left and 
right leads, according to the formula p n = (1^ — + 
ij:), where n = L,R. For R = (solid line) we have 
p L = p R for all biases. This curve is the same as seen in 
Fig. EKd). When R ^ these two polarizations depart 
from each other. The p L increases with R tending to 
reach the left lead polarization value pl = 0.23. In con- 
trast, p R assume negative values tending to —pr as R in- 
creases. This shows that even though we have a constant 
total current along the system [r L + 1£ = I R + 1^), its 
polarization can change across the system when R =/= 0, 
i.e., when there is a transverse magnetic field applied on 
the dot which coherently rotates the spin. 
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12 3 4 



parallel 




bias (meV) bias (meV) 

FIG. 3: (a)-(b) Spin-resolved currents 7f and and (c)-(d) 
current polarization for both parallel and antiparallel align- 
ments. In the P case we find 7f > 7j.. In contrast, for AP 
alignment 7f < 7j. in the bias window 1 — 3 meV (singly occu- 
pied regime) and 7j > when the bias exceeds the charging 
energy above 3 meV (doubly occupied regime). Panels (c)-(d) 
reveal a suppression and an enhancement of the current polar- 
ization in the singly occupied regime for the P and AP cases, 
respectively. Parameters: pl = 0.23, pr = 0.35, k B T = 200 
MeV, P = 10 fieV, R = 0, e = 0.5 meV and U = 1 meV. 

Figure [5] shows 1° against R. For R = 0, I L = l\ and 
= I R as expected. For R ^ these equalities disap- 
pear, with l\ and I R increasing and ij^ and I R decreasing 
with R. This leads to an enhancement of p L and \pn\ as 
seen in Fig. [4] Interestingly, there is a crossing point be- 
tween l\ and around R » To/4. So for this particular 
i? the total current becomes unpolarized in the emitter 
(left) lead and relatively high polarized in the collector 
(right) lead. This means that it is possible to change 
the current polarization from emitter to collector lead by 
preccssing the electron spin in the quantum dot. 

In the inset of Fig. [5] we show the total current against 
R. This curve resembles a typical Hanle resonance i 36 i 37 
Similarly to Ref. [3?]], here we can say that in the AP 
configuration and positive bias (i.e., with left being the 
emitter) the dot tends to be more up populated due to 
the majority up population in the emitter and the major- 
ity down population in the collector lead. On average a 
transverse magnetic field tends to increase the spin down 
component in the dot along the down magnetization of 
the collector lead. As a result the electron can more 
easily tunnel into the right ferromagnet and the current 
increases. 



C. Shot Noise 

Figure [5] shows the Fano factor, 7 = Sll/2sIl, against 
R in the AP configurations. The P alignment gives ap- 



FIG. 4: Current polarization against the bias voltage for dif- 
fering spin-flip rates in the AP alignment and in both left 
and right leads. For R = the polarization is the same on 
both sides. When R 7^ 0, though, the current polarization 
in the left side enhances, tending to the left magnetization 
degree pl- In contrast, the current polarization in the right 
side assumes negative values, tending to — pr. Parameters: 
p L = 0.23, pr = 0.35, k B T = 200 /zeV, Y = 10 /LteV, e = 0.5 
meV and U=l meV. 



antiparallel 

1.60-1 




1 1 ~2 3 ' 4 

R (10 2 meV) 



FIG. 5: Spin-resolved currents versus R. For R — we have 
I L = I R and 1^ = I R with i}. and /i being different from 
each other. For increasing R, l\ and I R grow while 7^ and 
I R decreases. For a particular 7? there is a crossing between I L 
and I L , thus resulting simultaneously in a non-polarized left- 
current and a polarized right-current. Parameters: pl — 0.23, 
p R = 0.35, k B T = 200 peV, To = 10 ^V, e = 0.5 meV, U=l 
meV and bias=4 meV. 

proximately insensitive Fano factor with respect to R. 
In the AP configuration, the Fano factor can be sup- 
pressed with R, reaching values below 0.5. This sup- 
pression can be further intensified by increasing the lead 
polarization parameters pl and pr. In particular, for 
fully spin-polarized leads (p — 1) AP aligned, the Fano 
factor reaches values close to 0.3 when double occupancy 
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is allowed (bias = 6 meV) and it attains 0.35 in the single 
occupancy regime (bias = 2 meV)^ For fully polarized 
leads in the P configuration the Fano factor remains at 
0.5 independently of R. 

A simple physical picture for this additional suppres- 
sion of 7 is as follows. Consider an up spin sitting on 
the dot. A second up spin trying to hop onto it is Pauli 
blocked, till the first electron tunnel to the collector lead 
or undergo a coherent spin-flip. If the spin-flip is fast 
enough (faster than the into/out tunneling processes) the 
first electron can return to its original state (via another 
spin- flip), instead of tunneling out of the dot. This blocks 
additionally the second up spin, consequently suppress- 
ing even further the noise. 

We note that for any R when we go from the single 
(bias = 2 meV) to the double (bias = 6 meV) occupa- 
tion regimes a reduction of the Fano factor is observed 
in both P and AP alignments (cf. solid black to solid 
gray lines and dashed black to dashed gray lines). This 
general feature was already predicted in Ref . [56| , where 
a diagrammatic formulation for the noise is derived. It 
is valid to mention that in the present study we have 
performed an S-matrix expansion (Appendix B) to the 
noise which could in principle be mapped into a Feyn- 
man diagrammatic formulation. Comparing our results 
with previous findings in the literature we observe a dif- 
ference between them in the single-occupancy regime^ 
Figure [7j shows a comparison for the shot noise obtained 



from Eq. (|3"0|) and from analytical results in Ref. [54 1 
(derived for pl = PR = and R = 0). While the second 
plateau [II in Fig. ([7])] coincides in both numerical and 
analytical cases, the first plateaus (I in the plot) do not 
coincided This disagreement is related to the Hartree- 
Fock factorization underlying our calculation. 

It is yet valid to note that without Coulomb interaction 
(U = 0) and for fully spin polarized leads antiparallel 
aligned, the Fano factor is given by 



_ 1 l /? 2 (5 
7 2 (1 + /3 2 ) 2 



/? 2 ) 



(34) 



where (3 = R/Tq. Eq. (|34j) is also found for a three tun- 
neling barriers junction^ Hence, for fully spin-polarized 
AP leads and non-vanishing spin-flip, the FM-QD-FM 
setup resembles a three-barrier geometry. 



0.9 
0.8 
0.7 
I 0.6 
g 0.5 

LL 

0.4 
0.3 



dashed line: parallel 
solid line: antiparallel 



black line: bias=2 meV 
gray line: bias=6 meV 




R (10 meV) 



FIG. 6: Fano factor against R for two fixed bias voltages cor- 
responding to single (bias = 2 meV) and double occupancy 
(bias = 6 meV) in both P and AP alignments. In the P case 
the Fano factor is approximatelly constante. In contrast, in 
the AP alignment the Fano factor can be suppressed due to 
spin-flip, reaching values below 0.5. In particular, for fully po- 
larized leads (pl = Pr = 1) the Fano factor can be strongly 
suppressed, assuming values close to 0.3 when double occu- 
pancy is allowed. Parameters: eo = 0.5 meV, (7=1 meV, 
F = 10 fieV and k B T = 200 fieV. 



We also observe that the spin-flip can change the current 
polarization when it flows from the emitter to the col- 
lector. It is even possible to have a polarized current in 
the collector while it is unpolarized in the emitter. We 
have derived an expression for the noise [Eq. (j3"0)) ] which 
exactly accounts for spin-flip but only approximately for 
Coulomb interaction. Finally, we found a suppression of 
the Fano factor to values well below 1/2 due to spin- flip. 

The authors acknowledge support for this work from 
the Brazilian agencies CNPq, CAPES, FAPESP and 
IBEM and the FiDiPro program of the Finnish Academy. 
The authors thank J. Konig and J. Martinek for helpful 
discussions. 



V. CONCLUSION 

Using the nonequilibrium Green- function technique we 
have studied the transport properties of a quantum dot 
coupled to two ferromagnetic leads. We consider both 
parallel (P) and antiparallel (AP) alignments of the lead 
polarizations. Coulomb interaction and coherent spin- 
flip are included in our model. We find that for distinct 
ferromagnetic leads the interplay between Coulomb in- 
teraction and spin accumulation translates into an en- 
hancement and a suppression of the current polarization 
p in AP and P cases, respectively, depending on the bias. 



APPENDIX A: DOT GREEN FUNCTIONS 

Here we present in some detail the calculation of the 
dot Green functions G r , G a , G > and G<, used in the 
current and noise expressions. The starting point is to de- 
rive an equation of motion for the contour-ordered Green 
function G cto -'(t,t') = — i(T c d CT (r)dJ.,(r')), and then, via 
analytical continuation rules, to determine these Green 
functions. After a straightforward calculation via equa- 
tion of motion we find 

G(r,r') = G°(r,r') 




equation for G 



1 2 3 

Bias voltage (meV) 



FIG. 7: Comparison of the shot noise for the present approach 
(numerical) and the one derived in Ref. [54| (analytical). 
While the second plateau (II) coincides in both approaches, 
the first plateau (I) differs. We attribute this contrast to 
the Hartree-Fock type approximation used in our calculation. 
Parameters: Pl ~ Pr = 0, ksT = 50 /ieV, To = 10 fj,eV , 
R = 0, e = 0.5 meV and U = 1 meV. 



dT 1 dr 2 G (T,T 1 )S(T 1 ,T 2 )G(T 2 ,r') 
+ J dnGV.rOtfGWfay), (Al) 

where the components of G' 2 ) and £ are 
GgW') = -t(T c nj(T)d t (T)d](r')), and S. y -(r,r') = 

Efa^M^H^W^r'), respectively. G° is 
the dot Green function without both the coupling to 
leads and the Coulomb interaction. 

Following the equation of motion expansion we find for 



G( 2 >(t 1; t') = G°( 2 )(r 1 ,r')+ J J dr 2 dr 3 G ^ (n,r 2 ) 
x S(r 2) T3)G(r 3 ,T , ) ) (A2) 
where G°( 2 )(ti,t') satisfies the identity 

(i^--e i -U)Gf\r,r') = S^d(t - - 

6 ij 6(t-lf)(£d i ).(Aa l ) 
In Eq. (|A2[) we have used the following approximations, — 
{Td-{t)d % {t)cl mi {t)d){t')) = 
(r4(t)d j (t)c fc ^(t)dt(t / )) = (dldi)(Tc kan (t)d}(t')) 
{Td\{t)dj{t)c kan {t)d){t')) = (dUi)(Tc kan (t)d}(t>)). 



Observe that Eq. (|A2[) closes the system of equations 
(|ATj) - (|A"2|) . Substituting ||X2]) into (|AT|) we find a Dyson 



G(r,r') = G°(r,r')+ / / dr^G^T, n) 

x S(r 1 ,r 2 )G(r 2 ,r'), (A4) 



where 



G°(T,r')-G (r,r')+ | dnG^r, n)[/G ( 2 ) (r 1; /)■ 

(A5) 

Applying analytical continuation rules in Eqs. (|A4|) and 
151) we find 



G r (t,i') = G Qr (t,t') + I I dtxdhG^^tx) 

x £''•(*!, i 2 )G r (i 2 ,i'), (A6) 

G 0r (t,t') = G 0r (t,t') + / «ftiG 0p (t > ti)^G°W(ti,* / ) ) 

(A7) 

and also the Keldysh equation 



G<(t,f) = / i dt 1 dt 2 G r {t 1 t 1 )^ < {tiM)G a (t 2l t'). 

(A8) 

Via Fourier transform of Eqs. (|A6|) and (|A7|) we obtain 



(A9) 



G r (e) = [G 0r ~ (e)-S r (e)]- 1 , 

G 0r (e) = G 0r (e) + G Qr (e)UG 0(2 \e), 
and of Eq. (|A"8|) we find 

G<( e ) = G r (e)S<( e )G a ( £ ), 

where 



G 0r (e) 



£-£l+iO+ 





2+t0+ 



(A10) 



(AH) 



(A12) 



and 



G 0(2)r (e) 



(n 2 ) (djdi) 

e-£l-f/+i0 + £-£l-C/+iO + 

-it. 



<«!> 



(A13) 



V £-£ 2 -C/+iO+ £-e 2 -C/+iO+ 

The retarded self energy is given by 



1 / r T + r x r T -r x 

2 I r T -r x r T + r x 



and the lesser self energy is defined as 

E< m = ^|(-i)(' +m )^(n i r^ + n ii r«). 



(A14) 



(A15) 
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APPENDIX B: S-MATRIX EXPANSION FOR 
THE NOISE 



According to Eq. (|21|) the noise is given in terms of the 

(2) 

four-operator Green functions g\ (r, t ) [i = 1,2,3,4]. 
To determine their equations of motion we develop an S- 

(2) 

matrix expansion as we illustrate below for g\ (r, r ). 
The first step is to transform the operators from the 
Heisenberg to the interaction picture, 

5l 2) (r,r') = l 2 (T c 5?L ) (r)d l (r)4 ff , t; ,(r')d J (r')), 

where the tilde denotes the operators in the interaction 
picture, i.e., 



di(t) = v\t,to)dt(t)v(t,t ), 



with 



,. , n rji i /' dt'H T (t') 

v(t,to) = Te Jt o 



(Bl) 



(B2) 



and a similar definition for the c operator. The oper- 
ator T is the time-ordering operator. The 5-matrix in 
g^\r, t') is defined as 



5 = yg-'/c dTtHriTi) 

Expanding S* we find 

n=0 



(B3) 



k 1 a 1 r) 1 i 1 



+ *fe 1 <, 1 ^ 1 (n)4 1 ^ I (n)] n ), 



(B4) 



where the lowest-order nonzero term in the expansion is 
that of n = 2. Since we assume nonintcracting leads, we 
can factorize the angle bracket in Eq. (|B4[) into a product 
of the lead and dot parts. We then apply Wick's theorem 
to the lead part. This results in 



9 ?>(t,t') = t* ka tl, a , J J d Tl dT 2 

E k-i^-K-i 



k 1 a 1 r] 1 k 2 a 2 r] 2 



ii 12=1,2 

( — *)(T , cCfe 2CT2 ^ 2 (T 2 )c^, (T , 1J ,(t')) 

(^(^^•(rO^Cn)^^)). (B5) 

In Eq. (|B5[) we have contracted C); icrim with c|. , and 

Cfc 2 (T2?j2 w ith cj. /cr , / . This is one choice among n(n — 1) 
possible contractions. Since all of them yields the same 



result, we simply multiply the chosen pairing by n(n— 1). 
This factor cancels part of the factorial n\ in Eq. (|B4p . 
thus resulting in the 5-matrix in the last angle bracket 
ofEq. ([B5]> . 

The first and second averages in (|B5[) give SkarjMcivi 
and Sk'cr'ri' ,k 2 ,o- 2 ,J72) respectively, so the sums over 
(&X, o"i, ?7i) an d (^2i f2) disappear. Defining 
9ka V (Ti,T) = -i(r c Cfe CT ,(ri)4 CT?? (r)) and gya'v'fa, r 2 ) = 

-i(7 7 c Cfc' -'77'( r ')Cfc' < r'77'( T 2)}) WG can rewrite (jB5j) as 



3*1 ,*2 — 1,2 



x fffan, (n , r)g k 'a' V ' fa , t' ) (T c dj (r)dj (r')dT ( n )dj a (r 2 ) ) • 

For the C/ = case the calculation is straightforward. 
By applying Wick's theorem in the four operators Green 
function, we find 

5 f ) (r ) r / )=«v / / dndn 



x Y,l(- i y iS " i (- i y 2S °' i 9^nfa,T) 9 k'^fa,T') 

iii 2 

x {Gu 1 (t, T\)Gji 2 (t , t 2 ) — Gi-j 2 

(B6) 

where Gu 1 {t, n) = — i(T c cij(r)dj (ti)), plus analogous 
definitions for the other Green functions. A similar cal- 
culation yields Eq. (f2"3"| for (t,t'). In the presence 
of the Coulomb interaction (U ^ 0) Eq. (|B6|1 is no 
longer exact and the full diagrammatic expansion should 
be considered in order to find an accurate noise expres- 
sion. However, this is a formidable task since it involves 
not only the usual many body expansion but also the 
analytical continuation of two and more particles Green 
functions. So as a first approximation we use Eq. (|B6[) 
even in the presence of the Coulomb interaction. 



APPENDIX C: "LINKED-CLUSTER THEOREM" 
FOR THE NOISE EXPANSION 



Equations (f22|) - (|25)) arc composed of what we call con- 
nected and disconnected terms. Here we show that the 
disconnected parts cancel identically the term 21^ in Eq. 
(f2"Tj) . Writing explicitly the disconnected term of Eq. (f2"2")) . 
we have 



9i 



?Lfar') = tlJl, a , J(-l)^(-l) i2Vi 



X y dTlG lll (T,T 1 )g ka11 (Tl,T + ) 

X J dT 2 G Tl2 {T / ,T 2 )gk'a'r,'fa,T' + ), (CI) 
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where the + sign on one of the r and t' is just re- On the other hand we can write the current as 
minder that the sequence of operators c k (T)di(T) and 



c !'erV ( T ')dj ( r ') m the main definition of 9^\t,t') (be- 
ginning of Sec. III-D) should be preserved during the 
following calculation. Applying the analytic continuation 
rules we obtain 



Fia,k-q(i{t, t)Fja' ,fe'j7'cr' (t , t'), (C2) 



with 



ii 

(t, ii ) 5fc % (*!,*) + Gt (t, h)g a kav (tut)} ,(C3) 

and a similar definition for i^v ,k'ri'a'{t' \t'). Similarly, 
from Eq. l|2"5)) , we have 



ii,i 2 =lj2 

X J dT 1 G lll (T,T l )g karl (T ll T + ) 

x y rfr 25fe v, I) -(r',r 2 )G l2 j(r 2 ,r /+ ). (C4) 
which, after analytic continuation, can be expressed as 

= -Cw 

F iaMa {t,t)F* a , tk , n , a ,{t',t')- (C5) 
sing the identities Eqs. (|24l) - ([25|) we obtain 

= -t ka t* k , a , 

Fl Ma (t,t)F ja ,^, a ,{}!j\ (C6) 

and 

F* aMa (t,t)F* a ,^, a ,(t',t>). (C7) 
From Eqs. and (J2TJ we note that 

kk'tjtj'ij 

{tkat k , a ,g^> c {t,t') 

-*fc<7*fc'<7'S2di^(*> *') 

-*fc<r*fe'<T'53difc(*) *') 

+ft.c. (C8) 



Using Eqs. JCT]), HC4J), JC6]) and jC7]) in Eq. flCS) we 
find 

({^(t),V(i')})disc = 2e 2 £ |t feCT | 2 |W| 2 

x [F^, k ^At',t')+FU',kwAt',t% (C9) 



(!„(*)) = 2eRcEl^| 2 ( F laMa (t,t). 



V2 



X [^^(t.tJ+K.fcvr (*,*)] (CIO) 

Squaring Eq. (|C10p and multiplying it by two, we find 
2(T v (t))(T v ,(t'))=2e 2 E |W| 2 |W| a 

x [^^(i'^O+^^v'^'O]- (en) 
Hence, Eq. (|C11[) cancels identically with (|C9[) . i.e., 

({T n (t), V(t')})disc - 2(/ r; (t))(V(t')) = 0. (C12) 



APPENDIX D: RECOVERING THE STANDARD 
FORMULA FOR THE NOISE 

To prove Eq. (p?Tj) we note that the Green functions 
appearing in Eq. (|30[) can be written as follow, 

T L G< = iT L G r {n L T L + n R T R )G a , (Dl) 



T h G> = -iT L G r [{l - n L )T L + (1 - n R )T H ]G a , (D2) 

r L (G r - g q ) = -ir L G r (r i + r fl )G a , (D3) 

Tr[r L G a r L G a + r i GT i G r ] = Tr[2r L GT L G a + 
r L (G r - G a )r L (G r - G a )].(D4) 

Now, defining the generalized transmission coefficients 



T LL = r L G r r L G Q 

T LR = r L G r T R G a , 



(D5) 
(D6) 



we can write the above set of equations [Eqs. (|E)l|) - (jE)4[> ] 
in terms of Tll and Tlr, 



r L G< =in L T LL +in R T LR , 



T L G> = -t(l - n L )T LL - i(l - n R )T LR , 



T L (G r - G r ) = -iT LL -iT 



LR, 



(D7) 



(D8) 



(D9) 



Tr[r L GT L G a + r i GT i G r ] = Tr[2T LL - 

(T LL + T LR )(T LL + Tlr)]- (D10) 
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Using Eqs. (|D7)) - (|DT01) in Eq. J30]) wc obtain 



Sll(u = 0) = 




nj,(l - n L )T LL + n L {l - n R )T LR 

+ (1 - n L )n L T LL + (1 - n L )n R T LR 

+ (n L T LL + n R T LR )[(l - n L )T LL + (1 - n R )T LR ] 

-n L (T LL + T LR )[(1 - n L )T LL + (1 - n R )T LR ] 

-(1 - n L )(T LL + T LR )[n L T LL + n R T LR ] 

-n L (l - n L )[2T LL - (T LL + T LR )(T LL + T LR )]}. 



The terms with Tll cancel out identically, and the above 
expression can be written as 



Sll(uj = 0) = 




x Tr{[?7. L (l - n L ) + n R (l - n R )]T LR (e) 
+ (n L ~ n R ) 2 T LR {e)[l - T LR (e)]}. (Dll) 

Denoting Tlr simply as T we arrive at Eq. (f5Tj) . 
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